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Abstract. We present a detailed analysis of the photometric properties of galaxies in the cluster A2163B at redshift z ~ 0.2. 
R-, I- and K-band structural parameters, (half light radius R c , mean surface brightness <yu> c within R c and Sersic index n) 
are derived for N ~ 60 galaxies, and are used to study their internal colour gradients. For the first time, we use the slopes of 
optical-NIR Kormendy relations to study colour gradients as a function of galaxy size, and we derive the Photometric Plane at 
z ~ 0.2 in the K band. Colour gradients are negligible at optical wavelengths, and are negative in the optical-NIR, amounting 
on average to -0.48 ± 0.06. Fhis result is in agreement with our previous measurements of colour gradients at intermediate 
redshifts, and imply a metallicity gradient in galaxies of ~ 0.2 dex per radial decade. Fhe analysis of the Kormendy relation 
suggests that its slope increases from the optical to the NIR, implying that colour gradients do not vary or even do become less 
steep in more massive galaxies. Such a result is not simply accomodated within a monolithic collapse scenario, while it can 
be well understood within a hierarchical merging framework. Finally, we derive the first NIR Photometric Plane at z ~ 0.2, 
accounting for both the correlations on the measurement uncertainties and the selection effects. The Photometric Plane at 
z ~ 0.2 is consistent with that at z ~ 0, with an intrinsic scatter significantly smaller than the Kormendy relation but larger than 
the Fundamental Plane. 
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<■£ 1. Introduction redshift up to z ~ 1, while its zeropoint changes accord- 

. . \ ing to what is exp ected for an old passively evolving stel- 

> Studies of galaxy populations in clusters at intermediate red- lar popu i ation (e.g. LStanford. E isenhardt and Dick insonlll998i 

shifts have been proved to be an effective tool to constrain | Ko dama et al.lfl998l here after KAB98). As shown by KAB98 

the mechanisms underlying their formation and evolution. A (see also iMerluzTet^Eool) . the evolution of the CM se- 

pure photometric approach involves the analysis of galaxy quen ce with redshift sets strong constraints on the origin of this 

Luminosity Functions (LFs), of galaxy colours, mainly via the re i a tion, resolving the well known age-metallicity degeneracy 

colour magnitude (CM) diagrams, and of the internal light dis- dWorfhev et al i ll 996) 
tribution of galaxies, i.e. their structural parameters. 

The galaxy LF and its dependence on the waveband 
carry information both on the stellar processes occurring 



The structural properties of galaxies at different redshifts 

provide a wealth of information by which galaxy evolution can 

, . . , be investigated, both through the study of the correlations and 

in galaxies and on the cosmological processes which deter- . . . . <- , 

. , . . _. rr — c — . " , — : — i , i through the analysis of the waveband dependence of such quan 

mine their mass function. Since Press & Schechter ( 1974), . . „. XTTT , .. . . „ . . 

the shape of the LF has been seen as a basic test of the- 
ories of structure formation and evolution. The luminosi- 
ties of early-type galaxies are also tightly correlated with 
their colours, through fhe colour magnitude (CM) rela- 
tion dVisvanathan & Sandagel 1 19771 iBower. Lucev and Ellisl The multi-waveband analysis gives information on the in- 



tities. The NIR light, in fact, is less sensitive to more recent star 
formation, to metallicity effects (though line blanketing) and to 
dust absorption, following more closely the luminous matter 
distribution. 

The mul 

Il992h . The slone of this relation shows little evolution with ternal colour gradients (CGRs) of galaxies, and therefore on 

the radial gradients of the properties of their stellar pop- 
Send offprint requests to: F. La Barbera ulations (SPs), such as age, metallicity and dust content. 

* Based on observations collected at European Southern Nearby galaxies have on average negat ive CGRs, their SPs 
Observatory (ESO 65.O-0251, 69.D-0653). Tab. 2 is fully available in becoming bluer toward the peripher y dPeletier et all Il990t 
electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr. IPeletier. Valentiin and Jame son. 1990). Due to the low evolu- 
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tion of colour gradients with redshift, metallicity seems to be 
the primary driver of the SP gradients in spheroida l galaxies, as 
demonstrated both by optical- optical studies ("e.g. lSaglia et alJ 
2000] iTamura and Ohtall2000|) and by opti cal-NIR CG mea- 
surements (La Barbera et all bo02l l2003al hereafter LBM02 
and LBM03a). Age gradient s, howeve r, are not fully ruled out 
by the present data dSaglia et all2 000. LBM03a), while the ef- 
fects of dust absorption still remain substantially unresolved. 
The presence of metallicity gradients can be well explained 
in the m onolithic formation scenario of early-type galaxies 
llLarsorifl974l) . due to the later beginning of a galactic wind 
in the inner galaxy regions. However, it can also be accommo- 
dated within the hierarchical merging framework, since metal- 
licity gradients are settled in disk galaxies at high redshifts, and 
more massiv e early-types form on average from the merging of 
larger disks jKauffmana 1 1 9961) . 

Structural parameters, such as the effective (half-light) 
radius R e , the mean surface brightness <fi> e within R e 
and the Sersic index n, are correlated by various relations, 
whose physical origins reside both in the properties of the 
SPs of galaxies and in the dynamical structure of these 
systems. For the early-type galaxies, one of these correla- 
tions is that between logR e and <fi> e , also known as the 
Kormendy relation (KR). As shown by Capaccioli, Caon 
and D'Onofrio (1992), early - type galaxies + bulges form 
two distinct families in the plane of effective parameters: 
that of the bright early-types, following the KR, and another 
'ordinary' family, whose properti es are more dispe rse and 
heterogeneous. Recently, Graha mand Guzminl ( 120031) showed 
that a real dichotomy does not exist between the two fam- 
ilies, but their observed properties are well explained by a 
systematic change in the profile shape with galaxy luminosity. 
The evolution of the zeropoint of the KR has been largely 
investigated at optical wavebands to constrain the major for- 
mation epoch of galaxy SPs a nd to perform the Tolman t est for 
the cosmological expansion ( S andage and Perelmuter , 1991; 
Pahre. Diorgovski and de Caryalhd 1 1 994 iBarger et al 1 
1998; Ziegl eretall Il999t ISandage and Lubinl [2001; 
Lubin and^Sandagei l200ll) ~"on the other hand, the slope 
of the KR is an interesting tool to gain insight in to the proper- 
ties o f the galaxy SPs as a funct ion of galaxy size JZiegler et alJ 
1 19991 lLa Barbera etaD 12003 hi hereafter LBM03b), which is 
a crucial predict ion of hierarchi cal merging scenarios. As re- 
cently shown bv lGrahaml (120021 hereafter GRA02), early-type 
galaxies also follow a three-dimensional relation between 
logR e , </U> e and logn, which is similar to the spectro- 
scopic Fundamental Plane (FP), once velocity dispersions 
are replaced by Sersic indices. This 'photometric plane' 
(hereafter PHP) has an observed scatter comparable to that of 
the Fundamental Plane (FP) relation, and carries interesting 
information on the physical properties of galaxies and their 
origin. 

In the framework of a project aimed at investigating the 
optical-NIR photometric properties of galaxies in the fields 
of ROSAT PSPC extended sources, we have obtained multi- 
waveband data (BVRIK) for a field of 5' x 5' centered at 
RA = 16hl5m48s and DEC = -06° 02' 10", at about 6.5' 
North (~ 1.3 Mpc at z ~ 0.2 with H = 70 Km s^Mpc" 1 , 



Q m = 0.3 and Q A = 0.7) 1 of the Ab ell cluster A 2163 at 
redshift z = 0.201 dArnaud et all Il994 . The coordinates of 
this field correspond to the center of a secondary extended 
emission in the X-ray map o f A2 163, known as A2163B 
(Elb az. Arnaud and Bohringeil 1 19951 hereafter EAB95), with 
an X-ray temperature KT > 6 keV. The cluster A 2163 (RA = 
16h 15m56s, DEC = -06°08'55") 2 is a very rich and com- 
plex structure of galaxies, which has been extensively studied 
for its exceptionall y hot X-ra y temperature, KT ~ 12 — 15keV 
("see iMarkevitch et alj fl996. and references therein), and for 
its huge radio halo , one of the most powerful known so far 
jFeretti et all EoOl. hereafter FFG01). Our data show the pres- 
ence of a significant excess of galaxies at redshift ~ 0.2 in the 
field of A2163B, with a NIR spatial distribution typical for a 
cluster of galaxies, with a main central overdensity and a sec- 
ondary structure at ~ 0.4 Mpc in the North-East direction. 

In the present work, we study the photometric and struc- 
tural properties of the galaxy population of A2163B, by dis- 
cussing the constraints implied by our results on different 
galaxy evolution scenarios. 

The layout of the paper is as follows. The data are presented 
in Sec. [2] while the luminosity density map, the colour magni- 
tude relations and the K-band LF of A2163B are studied in 
Sec. [5] Cluster members are then selected by the photometric 
redshift technique, as described in Sec. 0] The derivation of 
structural parameters is outlined in Sec. 5. Sec. 6 deals with 
the comparison of structural parameters and with the analy- 
sis of the internal colour gradients of galaxies. Sees. 7 and 8 
deal with the correlations between the structural parameters. In 
Sec. 7, we attempt for the first time to perform an optical/NIR 
comparative study of the log R e -< y. > e relation, while in Sec. 8 
we use the K-band data to analyze the PHP for the population 
of spheroids in A2163B. Discussion and conclusions follow 
in Sec. 9. Details on data reduction are given in Appendix A, 
while the estimate of the cluster redshift from the optical-NIR 
CM relations is detailed in Appendix B. The catalogue, in- 
cluding all the photometric properties of galaxies in the field 
of A2163B, that is BVRIK total magnitudes, colours, photo- 
metric redshifts and RIK structural parameters, is described in 
Section 2.3. 

2. The data 

2.1. Observations 

The data relevant for the present study were collected at the 
ESO New Technology Telescope and include BVRIK photom- 
etry centered on the field of galaxies A2163B. The data were 
obtained on August and September 2000 (hereafter run I), un- 
der partly non photometric conditions, and on April 2002 (run 
II), under photometric conditions. 

During three nights of run I, we obtained K-band imag- 
ing with the SOFI instrument (pixel scale 0.288"/pxl) for a 

1 We adopt this cosmology throughout the paper. With these pa- 
rameters the age of the universe is ~ 13.5 Gyr, and the redshift 0.2 
corresponds to a lookback time of ~ 2.5 Gyr. 

2 The coordinates refer to J2000 and define the peak of the X-ray 
emission (EAB95). 



F. La Barbera et al.: The cluster of galaxies A 2163B. 



3 



Table 1. Observations. The notes refer to the sky conditions: 
NP and P denote non-photometric and photometric conditions, 
respectively. 
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field of 5' x 5'. The sky conditions were photometric for two 
nights and non photometric for the third night. A total of 120 
exposures were collected, with a DIT of 6 s, NDIT =10, and 
a dithering box of 18". Due to the non photometric conditions, 
only 85 exposures were retained, re sulting in a total integ ration 
time of 5100 s. Standard stars from Persson et al. (1998) were 
observed during each night, each at five different positions on 
the chip. In order to perform the illumination correction (see 
Sec. lA.2> . we also observed a standard star at different positions 
on a 5 x 5 grid across the frame. In the same period, we also 
obtained V-, R-, and I-band photometry with the EMMI instru- 
ment for a region of 8.6' x 8.6' (pixel scale 0.267"/pxl), under 
non photometric conditions. On April 2002, the field A2163B 
was re-observed in B, V, R and I band with EMMI under photo- 
metric conditions. The relevant information on the observations 
are summarized in Tab. ^ The total exposure times amount to 
600, 3000, 1200 and 1200 s for the B, V, R and I bands, respec- 
tively. 



2.2. Analysis 

The data were reduced by using FORTRAN routines devel- 
oped by the authors and IRAF packages 3 . Details are given 
in Appendix IA.1I and IA.2I for the optical and the NIR im- 
ages, respectively. Photomet ry was derive d by using the soft- 
ware S -Extractor dBertin & Arnoutl| 19961) . For the optical im- 
ages, S-Extractor was run in the so-called double-image mode, 
using the R-band image for the detection of the sources, and 
the other images only to measure their fluxes. Particular care 
was taken so that closed objects were deblended in the same 
way both in the optical and in the NIR images. The BVRI 
and the K-band catalogues were cleaned by spurious detections 
and cross-correlated. For each object, we computed colour in- 
dices by using magnitudes within a fixed aperture of diame- 



ter 3.5", corresponding to ~ 12 kpc at redshift z = 0.2. Since 
the images have different seeing, the aperture magnitudes were 
corrected to th e same FWFJM of 1.2", following the proce- 
dure detailed in La Barbe ra et al . (2003c). The aperture correc- 
tions amount to -0.04, 0., 0.02, 0.02, and 0.04 mag for the 
B, V, R, I and K bands, respectively. Total magnitudes were 
obtained with the S-Extractor mag-auto parameter, which pro- 
vides magnitud es in an ad aptive aperture of diameter a ■ 
where tk is the iRronl fl998) radius, and a = 2.5, for which 
the Kron magnitude encloses 94% of the object source. Total 
magnitudes 4 were then obtained by subtracting 0.06 mag from 
the Kron magnitudes. The uncertainties on galaxy magnitudes 
were obtained by adding in quadrature the errors on the photo- 
metric zeropoints to the uncertainties estimated by S-Extractor. 
Galaxy/star separation was performed on the basis of the stel- 
lar index SG of S-Extractor, defining as stars the objects with 
SG > 0.97. The completeness of the K-band photometry, which 
is used in Sec. 13.21 fo r the LF analysis, was derived by fol - 
lowing the method oflGarilli. Maccagni and Andreonl jl999). 
This method consists in the determination of the magnitude at 
which objects start to be lost since they are below the bright- 
ness treshold in the detection cell. To estimate the complete- 
ness limit, the magnitudes in the detection cell were computed 
and compared to the total magnitudes. This approach allows 
to define a conservative value for the completeness magnitude, 
corresponding to a completeness percentage of ~ 100%. The 
completeness magnitude in the K band is Kj = 19. Down 
to this limit, the BVRIK catalogue includes N = 317 galax- 
ies. The completeness of the BVRIK catalogue with respect 
to the K-band catalogue is ~ 98% at K T = 17.5, falling to 
~ 90% at K T = 19. At K T = 17.5, the typical photomet- 
ric accuracies on the B-, V-, R- and I- band data amount to 
0.38,0.17, 0.11 and 0.1 mag, respectively. For the optical cata- 
logues, we also computed the fraction of sources (fs) for which 
S-Extractor provides reliable magnitude estimates 5 as a func- 
tion of the K-band magnitude. For the B band, fs amounts to 
~ 90% at K T = 17.5, dropping to ~ 0.5 at K T = 19. For the V, 
R and I bands, we have f s (K T = 17.5) = 0.95,0.98,0.98 and 
f s (K T = 19) = 0.76, 0.84, 0.88, respectively. 



2.3. The catalogue 

The BVRIK photometric catalogue for the N = 131 galaxies 
with photometric redshift estimate (see Sec. 4) 6 is shown in 
Tab. 2. 

The catalogue is organized as follows. Two rows corre- 
spond to each galaxy. The first one provides the measured 
quantities, while the second one gives the relative uncertainties 
(one sigma standard intervals). All the photometric quantities 
have been corrected for galactic extinction by using the 
photometric redshift estimates and adopting E(B - V) = 0.415, 
with the extinction law from SFD98. The typical uncertainties 
on the structural parameters are discussed in Sec. 16. II 



3 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



4 In the following, total magnitudes will be indicated by a pedex T. 

5 We adopted an upper limit of 100% for the uncertainty on the flux. 

6 The full catalogue will be available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr. 
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Fig. 1. K-band image (left panel) and luminosity-weighted map (right panel) for the field of galaxies A2163B. The image is 
centered at RA = 16h 15m 48s and DEC = -06° 02' 10", and was smoothed for the plot by a gaussian convolution. The scale 
and the orientation of the image are shown in the lower left corner of the left panel. In the same panel, the crosses indicates the 
positions of the radio sources in the field. In the right plot, the circles enclose the regions that were used to calculate the LF in 
the higher density environment (see Sec. l3.2t . 



Column 1 : running number of the catalogue. 
Column 2: photometric redshift. 

Column 3, 4: right ascension and declination referred to J2000. 
The astrometric solution was computed by using a list of stars 
from the USNO catalogue. The rms of the residuals to the 
astrometric solution is 0.2". 

Columns 5-9: B-, V-, R-, I- and K-band total magnitudes. 
Columns 10-13: optical-NIR galaxy colours in the aperture of 
3.5" corrected for seeing effects. 

Columns 14-16: mean surface brightness <fi> e within the 
effective radius R e , logarithm of R e (in arcsec) and Sersic 
index n for the R-band. 

Columns 17-19: the same of columns 14-16 for the I-band. 
Columns 20-22: the same of columns 14-16 for the K-band. 
Columns 23, 24: axis ratio and position angle in degree 
obtained from the K-band surface photometry. 



3. The cluster A 2163B 

Fig-Hplots the K-band image of A 2163B and the luminosity- 
weighted density map for the N = 105 galaxies brighter than 
Kj = 17. The figure clearly shows the presence of a cluster of 
galaxies, with a main over-density in the center and a secondary 
structure at a distance of ~ 2' (~ 0.4 Mpc) from the main peak. 
We note that, on the basis of field galaxy counts (Sec. l3.2t . only 
~ 15 galaxies out of 105 are expected to be field contaminants 
in the area of the K-band image, confirming that the features in 
Fig- Upright panel) are highly significant. It is very interesting 
to note that both the X-ray image of A 2163B and the 20 cm 
radio map show a spatial structure very similar to the luminos- 



ity density image (see Fig. 7 of FFG01). In Fig^ we also show 
the positions of the two radio sources found by FFG01. One of 
this turns out to coincide with an early-type 7 galaxy (possible 
merging) in the main clump, while the second one is very close 
to a bright early-type in the second overdensity region. 

In order to derive the redshift z c of the galaxy popula- 
tion of A2163B we compared the BVRIK photometry with 
predic tions of stellar population models from the GISSEL00 
code (IBruzual & Charlotl Il993l) . In particular, we used (1) 
the optical-NIR colour magnitude (hereafter CM) relations 
(Sec. ED, (2) the K-band LF (Sec.E3 and (3) the photomet- 
ric redshift technique (Sec.0}. We point out that, although it is 
known that broad-band colours are degenerate in age, metallic- 
ity and redshift, and that uncertainties exist among stellar pop- 
ulation models (e.g. ICharlot. Worthev and Bressanl ll996). the 
large wavelength baseline provided by our photometry (from 
U to NIR restframe) and the use of the GISSEL synthesis code 
allow to obtain a reliable redshift estimate of galaxies at in- 
termediate redshift by using pur ely photometric techniques, 
as shown by other studies (e.g. |de Propris & Pritchel [l 998 , 
Busarello et al. 2002[ lLa Barbera et al.l2003ci) . Another impor- 
tant point for estimating z c is the effect of the galactic redden- 
ing E(B - V). The main sources for the reddening of the Milky 
Way a re based on HI measurements from Burstein & Heiles 
dl984 hereaft er BH84) and on COBE/DIRBE and IR AS/ISSA 
FIR data from lSchlegel. Finkbein er and Davisl ill 9981 hereafter 
SFD98). GeneraUy, the estimates of E(B - V) from SFD98 are 
systematically higher than those of BH84. For A 2163B, BH84 
and SFD98 predict very different extinction values, E(B - V) = 
0.18 and E(B - V) = 0.415, respectively. In the next section, 

7 Sersic index greater than 2. 
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Fig. 2. Optical-NIR CM diagrams for galaxies in the field of A 2 1 63B . For each panel, the slope and zeropoint of the CM relations 
are shown in the lower left corner, while the error bars from the left to the right denote the averaged uncertainties on colours and 
magnitudes at Kj = 15 and Kt = 18, respectively. 



we obtain a simultaneous estimate of the galactic reddening 
as well as of the cluster redshift from the optic al-NIR colour 
magnitude relations. This procedure gives an extiction value 
which is in remarkable agreement with that of SFD98 and a 
redshift estimate of z c ~ 0.2, which concides within the uncer- 
tainties with that obtained from the K-band LF (Sec. l3.2t . for 
which extinction effects are negligible. The difference between 
BH84 and SFD98 reddening values, therefore, does not affect 
the present analysis. 

3.1. Colour magnitude relations 

Fig.|2shows the optical-NIR CM diagrams, not corrected for 
galactic extinction, for all the galaxies brighter than Kt = 18.6 
in the field of A2163B. Each diagram shows a sharp sequence 
of red galaxies, following the CM relation: 

N-K = F Kr +A-(K T -K r ), (1) 

where N is the optical waveband, K r is a reference K-band 
magnitude, A is the slope and Tk, is the zeropoint of the re- 
lation, that is the N - K colour at Kt = K r . The values of Fk, 
and A, that are shown in Fig.|2j were estimated by applying a 
least square fit, minimizing the rms of the residuals in colour, 
cr N _ K , to the relation. In order to reduce the effect of outliers 



and field contaminants, the value of <tn~k was co mputed b y us - 
ing the bi-weight statistics llBeers. Flvnn and GebhardA ll990). 
Moreover, for comparison with previous works (KAB98), the 
fits were performed by considering only the galaxies within the 
three brightest magnitudes of the red sequence, corresponding 
to K T < 16.5. We note that, down to this limit, the BVRIK 
photometry is complete (see Sec. l2.2> . and that, on the basis of 
field galaxy counts (see Sec. l3.2> . only ~ 10% of the galaxies 
are expected to be field contaminants. Despite the large uncer- 
tainties, the values of A show that the slope of the CM rela- 
tion increases systematically from about -0.1 for the B - K 
sequence to about -0.04 for the I - K relation. The slopes of 
the V - K, R - K and I - K sequences can be straightly com- 
pared with those predicted from a pure mass-metallicity model 
by KAB98 at the redshift of A2163B (z ~ 0.2, see below): 
A v -K = -0.075, A R _ K = -0.062, and A^ K = -0.05 (see their 
Fig. 4). We note, in fact, that the CM slopes of KAB98 were 
normalized to a physical aperture of ~ 10 kpc, which is very 
similar to that adopted in the present work, and, therefore, the 
effects of galaxy colour gradients do not affect the compari- 
son of CM slopes. Our results fully agree with the values of 
KAB98. 

The zeropoints of the CM relation were used to obtain a si- 
multaneous estimate of the cluster redshift z c and of the galac- 
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Fig. 3. K-band LF of A2163B. The two curves indicate the 
best-fit Schecter functions obtained by fitting both a and K* 
(solid line) and by fitting only K* with a fixed to -0.9 (dashed 
line). In the last case, the three brightest cluster galaxies were 
excluded from the fit: the first bin was not considered while 
the galaxy counts for the second bin are marked by the empty 
circle. Contours in the smaller panel denote Icr ,2cr and 3<x 
probability levels, respectively. 

tic reddening in the direction of A 2163, fitting the values of 
T Kr with the colours expected from old stellar population mod- 
els (see Appendix B for details). In this approach, we took 
into account the degeneracy among metallicity and the other 
stellar population parameters exploiting the fact that the CM 
relation mainly ori gins from a mass-metallicity relation (see 
Mer luzzi et al.l Eo03l and references therein). This procedure 
gives z c = 0.215 ± 0.015 and E(B - V) = 0.41 + 0.02. It is 
quite remarkable that (1) the estimate of E(B - V) is in very 
good agreement with that of SFD98, E(B - V) = 0.415, and (2) 
the value of z c coincides with that obtained from the K-band 
LF, which is virtually unaffected by galactic reddening. For this 
reason, we adopted in the following E(B — V) = 0.415, and ap- 
plied the corresponding absorption corrections to our photom- 
etry 8 . Repeating the zeropoint fitting with E(B - V) = 0.415 
gives z c ~ 0.21, which is fully consistent with the redshift of 
the main cluster structure (A 2163 ) at z ~ 0.2. 

3.2. Luminosity function 

The evolution of the characteristic magnitude K* of the K- 
band LF of cluster galaxies at z < 1 is consistent with the 
behaviour of an old passive stellar population with formation 
redshift Zf e [2, 3] and solar metallicity. This result has been 

8 For E(B - V) = 0.415, at z ~ 0.2, and practically for any galaxy 
spectral model, we have A B = 1.630, A v = 1.281, A R = 1.049, Ai = 
0.788 and A K = 0.152. 



obtained bv ldePropris et alJ i 19991 hereafter dPS99) fitting a 
Schechter function (hereafter SF) with fixed slope a = -0.9 
to the K-band counts of cluster galaxies with K T - K* < 3, 
and turns out to be independent of cluster richness or X-ray 
luminosity. 

We used these results in order to obtain a further estimate 
of the redshift of A2163B. The K-band LF was obtained by 
correcting our number counts with field galaxy counts from 
the C alar Alto Deep Imaging Survey (CADIS, see lHuang et alJ 
2001 ), which is the largest medium deep K-band survey to date, 
with a total area of 0.2 deg 2 and a completeness magnitude of 
19.75 mag. Since our K-band photometry is complete down 
to K T = 19 (see Sec. 12.211. the CADIS data are the most suit- 
able dataset for the estimate of field galaxy counts. The LF of 
A2163B is shown in Fig. [3] The error bars take into account 
Poissonian uncertainties on both field and cluster counts, while 
the different curves show the SF fits of the LF. In order to ac- 
count for the finite size of the magnitude bins, the fits were per- 
formed by the convolution of a SF with the bin size. Following 
dPS99 and lde Propris & PritcheJ dl998l hereafter dPP98), we 
fitted the galaxy counts down to Kt = 18. (~ K* + 3, see be- 
low) by a SF with a = -0.9, excluding the three brightest clus- 
ter galaxies. The fit provided a good description of the counts 
of bright galaxies 9 , giving the following estimate of the char- 
acteristic magnitude: Kg 9 = 15.3 ± 0.20. In order to estimate 
z c , we minimized the function: 

* 2 = [K5. 9 (z)-KS. 9 ] 2 , (2) 

where K* g (z) is the value of K* expected at redshift z, and was 
obtained by adding to the value of the characteristic magnitude 
of the Coma cluster (10.9 + 0.2, see dPP98) the luminosity dis- 
tance term and the K+E corrections, derived from GISSEL00 
spectral models with a Scalo IMF and an exponential SFR. 
It turned out that, by using models with different metallicities 
(Z/Z = 0.5, l.,2.) and formation redshifts (Zf > 1.5), the es- 
timate of z c does not vary significantly: z c = 0.20 + 0.02. This 
result is in very good agreement with that obtained from the 
CM relations in Sec. 13.11 and, therefore, still shows that the 
redshift of A2163B is consistent with that of the main cluster 
structure A 2163. We found that the estimate of z c is very robust 
with respect to the field subtraction procedure and to the value 
adopted for the galactic reddening, due to the low sensitivity 
of the K-band photometry to dust absorption. For example, by 
using E(B — V) = 0, we obtain z c = 0.225, and neglecting field 
counts gives z c = 0.21. 

In order to study the faint end of the LF, we performed 
a SF fit for Kr < 19 by treating both a and K* as free pa- 
rameters. As shown in Fig. [3] the SF, whose best-fit parame- 
ters are a = -1.22 ± 0.15 and K* = 14.45 + 0.50, provides 
a good description of the LF of A2163B. The cluster counts 
show some deviation from the SF fit only for the two bins at 
Kt = 15.75 and Kt = 17.25. Although this could be an in- 
dication of a bimodal behaviou r of th e LF, as found by previ- 
ous studies (see iMercurio et aljE003l and references therein), 

9 The only two bins for which galaxy counts deviate by more than 
l<x from the model are at Kt > 17. We note that limiting the fit range 
to K T < 17 does not affect our results, giving K^ 9 = 15.34 ± 0.2. 
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Fig. 4. K-band LF of A2163B in the low- and high- density 
environments. The points and the error bars corresponding to 
the high density region have been shifted by +0.1 mag for sake 
of semplicity. The curves show the best fitting Schechter func- 
tions. 

the deviations are only at the levels of 1.6<x and 1.7 <x, respec- 
tively, and are therefore not very significant. The value of a 
is in very good agreement with the value of -1.18 found by 
lAndreonI J200ll) for a rich cluster of galaxies at z = 0.3 in 
an area of ~ 1.3 Mpc 2 , which is similar to that analyzed in 
the present study (~ 1 Mpc 2 ), and is also consistent with val- 
ues of the faint end slope of the NIR LF for the Coma cluster , 
a = -1.0 (dPP98) and a = -1. 3 jAndreon & Pellol l2000l) . 
and for field galaxies, a = -1.09 jKochaneTetalTBoOll) . We 
also analyzed the dependence of the LF of A2163B on envi- 
ronment, by performing the SF fits separately for the galaxies 
inside the higher density region, defined by the circles in Fig.^ 
and for the galaxies outside that region. The two LFs are shown 
in Fig.|4]with the corresponding SF fits. Interestingly, we found 
that the values of a vary from -0.96 + 0.18 (higher density) to 
-1.42 + 0.2 (lower density), implying a steepening of the faint 
end slope at lower densities, in agreement with the findings 
of other studies (see Haines et al. 2003, submitted, and refer- 
ences therein). This steepening is due to the fact that, as shown 
in Fig. |4] the lower density environment is characterized by a 
higher fraction of faint galaxies, with a significant difference 
(> 2cr) in the low- and high- density counts at Kj > 18. 

4. Photometric redshifts 

In order to select cluster members, we used the photometric 
redshift technique, which has been proven to be an effective 
tool to estimate the redshift of cluster galaxies at intermedi- 
ate redshifts (e.g. lBusarello et alfeOollLa Barbera et alJ2003"cl 
and references therein). Photometric redshifts were estimated 
according to the Spectral Energy Distribution fitting method 



z =0.198±0.024 

C 
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Fig. 5. Distribution of photometric redshifts. The solid line and 
the dashed region indicate the estimate of the cluster redshift, 
z c , and the relative uncertainty (l<x standard interval). Dashed 
lines enclose the range of z p , Az, used to select cluster mem- 
bers. 



(see IMassarotti. lovino and Buzzonil l200lt IMassarotti et al.1 
12001 1 and references therein). In order to achieve a reason- 
able accuracy, we considered galaxies with signal-to-noise ra- 
tio S/N > 5 in at least three bands, limiting the sample to the 
N = 131 galaxies brighter than Kj = 17.5. Since the percent- 
age of galaxies 10 at redshift z > 1 with K < 17.5 is expected 
to be negligible (< 2%), we looked for redshifts in the range 
Z € [0.0, 1.0] with a step of 0.01, imposing that at a given red- 
shift galaxy templates were younger than the age of the uni- 
verse in the adopted cosmology. We used the GISSEL00 code 
in order to produce galaxy templates with a Scalo IMF and an 
exponential SFR e _,/T . The colours of E/S0, Sa/Sb, and Sc/Sd 
spectra were modeled by choosing r = 1,4 and 15 Gyr, re- 
spectively, while early-type galaxies with different metallici- 
ties were described by using E/S0 models with Z/Z =O.2, 0.4, 
1 and 2.5. The differential dust extinction of the Milky Way was 
included in the comput ation of model colours by adopti ng the 
extinction curve of Cardelli, Clavton a nd Mathisj (I198SJI) and a 
colour excess of E(B - V) = 0.415 (see Sec. 13 . II . The un- 
certainty on the photometric redshift 6z was estimated by per- 
forming numerical simulations, taking into account the mea- 
surement errors on galaxy colours. The mean value of 6z is 
~ 0.08, varying from ~ 0.05 at K T = 15 to ~ 0.1 at K T = 17. 

As shown in Fig. |5] the distribution of photometric red- 
shifts is dominated by the peak around z ~ 0.2, indicating that 
most galaxies with Kj < 17.5 are actually cluster members. 
Applying the bi-weight statistics, we obtained the following es- 
timate of the cluster redshift: z c = 0.198 ± 0.024, in agreement 
with the values obtained in Sec l3.ll and l3.2l Galaxies with pho- 



10 This estimate was obtained by using a Pure Luminosity Evolution 
model (see e.g. Pozze ttiTBruzual and Zamoranill993) . 
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tometric redshift in a range of 25z around z c were defined as 
cluster members, resulting in a final list of N = 102 objects 11 . 
This value is in good agreement with that predicted from the 
LF of cluster galaxies, N = 1 12 (see Sec. l3.2> . 

5. Surface photometry 

Surface photometry was derived for all the galaxies brighter 
than Rt = 19.5 in the R and I bands, and for the galaxies 
brighter than Kj = 16.5 in the K band. These selection cri- 
teria were chosen in order to obtain reliable structural param- 
eters from the A2163B photometry, and were the results of 
numerical simulations performed as described in LBM02. A 
further selection was performed a posteriori by excluding few 
faint galaxies in the I and K bands because of the strong con- 
tamination from closed stars. The resulting samples include 
N = 64, 62 and 58 galaxies for the R, I and K bands, respec- 
tively, with N = 62 objects in common between the R and I 
bands, and N = 58 galaxies in common between the R and K 
bands. Structural parameters for these galaxies are given in the 
catalogue (see Sec. 2.3). 

Galaxy images were fitted with 2D models convolved with 
the PSF, as detailed in LBM02. As shown in that paper, the 2D 
approach allows reliable structural parameters to be obtained 
for galaxies at intermediate redshift (z ~ 0.3) from ground 
based data taken under ordinary observing conditions (seeing 
FWHM ~ 1.0" and pixel scale ~ 0.3"/pxl). Galaxy models 
were parametrized by the Sersic law: 

I(r) = Io-exp(-b-(r/R e ) 1/n ), (3) 

where r is the equivalent radius, R e is the effective (half- 
light) radius, Io is the central surface brightness, n is 
the Sersic index, and b is a constant (b ~ 2n - 1/3, 
see ICaon. Caoaccioli and D'Onofrioll^gsl) . For each galaxy, 
nearby objects were masked interactively, while overlapping 
galaxies were fitted simultaneously. The x 2 minimization was 
performed by the Levenberg-Marquardt method, deriving a to- 
tal of six best-fit parameters: the center coordinates, the effec- 
tive radius, the central surface brightness, the position angle PA 
and the axis ratio b/a. The mean surface brightness within R e , 
</i> e , and the total magnitude mr were computed from the 
fitting paramet ers by using the well k nown properties of the 
Sersic law (see lCiotti and Bertinl l999). The typical uncertain- 
ties on galaxy parameters were estimated by comparing struc- 
tural parameters between different wavebands (see next sec- 
tion). The PSF models were obtained by a multi-Gaussian ex- 
pansion of the images of bright unsaturated stars in the cluster 
field. For the R and I bands, all the stars were fitted simulta- 
neously, resulting in a single PSF model. For the K band, the 
PSF turned out to vary across the field, showing significant de- 
viations from the circular shape in the North-East corner of the 

1 1 We note that in the analysis of the optical-NIR colour gradients we 
consider only the N ~ 60 galaxies brighter than K T = 16.5. Preceding 
as described in lLa Barbera et alll2003cf) . we estimated that at this lim- 
iting magnitude, in the redshift range adopted to select cluster mem- 
bers, only ~ 1% of the galaxies are expected to be field contaminants, 
proving the reliability of the cluster members selection for the surface 
photometry analysis. 



frame. Since star images were available at different positions 
across the field, with an almost uniform coverage, we adopted 
for each galaxy the PSF model obtained from the closest star. 
PSF distortions were treated as in LBM03a, by adopting gaus- 
sian functions with elliptical isophotes. 

Fig. [6] plots, as an example, the K-band 2D fitting of some 
galaxies in the cluster field. The models turned out to give a 
very good description of the galaxy images. 

6. Optical-NIR internal structure of cluster 
galaxies 

The waveband dependence of the structural parameters car- 
ries information on the differential properties of SPs inside 
galaxies. This subject was investigated in LBM03a, who first 
attempted to analyze the evolution of the UV-NIR structural 
properties of both the populations of disk dominated (Sersic 
index n smaller than 2) and spheroidal galaxies up to redshift 
z ~ 0.6, by using a large sample of cluster galaxies (N = 270). 
Following that work, we studied the internal optical-NIR struc- 
ture of galaxies for A2163B by (1) a straight comparison of 
structural parameters (Sec. 16. 1> and (2) estimating the internal 
colour gradients of galaxies (Sec. 16. 2> . Since the number of 
disks in common between the R and K bands turned out to be 
negligible (N = 5 out of 58), we did not attempt 12 to study the 
waveband dependence of their structural properties. The fol- 
lowing results will refer, therefore, to the spheroidal population 
ofA2163B. 

6.1. Comparison of structural parameters 

Fig. shows the comparison of structural parameters for the 
N = 62 galaxies in common between the R and I bands. 
Differences are always computed by subtracting the quanti- 
ties at longer wavelengths from those at shorter wavelengths. 
The following quantities are considered: effective radii and 
mean surface brightnesses, Sersic indices, and the combina- 
tion of R e and </-i> e which enters the Fundamental Plane: 
kpp = log R e -0.3- <yu> e - We note that the differences of mean 
surface brightnesses, A<fi > e , were obtained by subtracting the 
R-Igalaxy colours, thatis: A</i> e =</i>g - <n>\ -(R-I). 
The mean value, m, and the standard deviation, <x, of each dis- 
tribution are shown in the figure, while in the lower right panel 
we also plot the correlation between AlogR e and A < jj > e . 
Optical parameters are in remarkable agreement, with mean 
differences that are fully consistent with zero for each quan- 
tity. This is consistent with the fact that the R and I bands at 
z ~ 0.2 correspond approximately to V and R restframe, sam- 
pling, therefore, a very similar spectral region for early-type 
galaxies. For this reason, we used the distributions of Fig.fJJto 
estimate the uncertainties on log R e , <yu> e and log n, by com- 
puting the covariance matrix of the differences between the R- 
and I-band parameters. The standard deviations of the distribu- 
tions amount to ~ 25% in R e , ~ 0.45 mag/arcsec 2 for </i> e 

12 Disk galaxies are included in the analysis of Sec. 16. II and 16.21 
However, the results are unchanged by excluding these galaxies from 
the samples. 
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Fig. 6. 2D fits of galaxies in the K-band image. For each plot, from left to right, the galaxy image, the 2D model and the residual 
map are shown. The size of the horizontal white bar is 5". 

and ~ 18% for n. We note, however, that due to the well known 
tight correlation between th e uncertainties on the effectiv e pa- 
rameters (see e.g. Jeirgense n. Franx an d Kiaereaard 1995), the 
kpp values have a much smaller dispersion, only 6%. 

Differences between R- and K-band structural parameters 
are shown in Fig. [H] We note that the differences of mean sur- 
face brightnesses, A < ft > e , were obtained by subtracting the 
R - K galaxy colours, that is: A < fi > e =< H >f - < fi >f 
-(R - K). The mean values of the distributions show that clus- 
ter galaxies are more concentrated in the NIR than in the op- 
tical, having NIR effective radii which are on average ~ 40% 
smaller than in the R-band (AlogR e = -0.16), and Sersic in- 
dices slightly larger in the K band, although the last difference 
is not very significant. The mean difference of A < fi > e can 
be fully explained by using the definition of total magnitude 
mx = -2.5 log(27r) - 5 log R e + < fi > e , and computing the dif- 
ference between R- and K-band total magnitudes: 

Rx - Kt = -5AlogR e + <fi>f - <fi>f. (4) 

Since A <fi> e ^<fi>f - <fi>f -(R T - K T ), we obtain the re- 
lation A <n> e = 5A log R e , which, for AlogR e = -0.16, gives 
A <fi> e ~ 0.8, in very good agreement with the mean value 
shown in Fig. [8] The same procedure allows the mean differ- 
ence of the kpp variable to be explained. 
The dispersion of the distributions in Fig. [H] is due (1) to the 



uncertainties on both the R- and K-band parameters, and (2) to 
the intrinsic scatter of the optical-NIR properties of galaxies. 
We verified by numerical simulations (performed as described 
in LBM02) that the uncertainties on the K-band structural pa- 
rameters are as large as or smaller than those estimated for the 
optical wavebands. This suggests that the dispersions of the his- 
tograms in Fig. [3] have mainly an intrinsic origin. In Sees. Q 
and [8] in order to describe the uncertainties on logR e , <fi> e 
and logn, we adopted the same covariance matrix for both 
the optical and the NIR data. By looking at Fig. [8] we stress 
again the small dispersion in the kpp variable, that amounts to 
~ 0.08 dex. 

6.2. Colour gradients 

The internal colour gradient of galaxies can be computed from 
the values of the effective radii and of the Sersic indices in the 
different wavebands, as described in LBM02 (see their Eq. 4). 

The distributions of colour gradients for the galaxies of 
A2163B are shown in Fig.|5J we wrote the R - I and R - K 
colour gradients as V-R and V-Krestframe gradients, respec- 
tively. By using different galaxy templates from the GISSEL00 
code, we verified, in fact, that the conversion from R - 1 (R - K) 
at z = 0.2 into V - R (V-K) at z = is independent of the spec- 
tral type (for Zf > 1), with lcr variations of 0.01 and 0.03 mag 
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Fig. 7. Comparison of structural parameter in the optical. For each panels, differences are computed between R and I bands. The 
correlation between effective parameters is shown inside the lower right panel. 



between different spectral models for the optical-optical and 
the optical-NIR colours, respectively. As expected on the ba- 
sis of the structural parameters, galaxies show optical-optical 
colour gradients which are fully consistent with zero, while the 
V-K colour gradients are negative for most galaxies, implying 
that their SPs are on average redder in the center. The mean 
value of the optical-NIR gradient is -0.48, which implies that 
galaxies become ~ 0.5 mag bluer per decade of radius toward 
the periphery. The value of grad(V - K) is fully consistent with 
the colour gradient estimate of LBM03a (~ - 0.4 mag/dex). 

7. Optical-NIR Kormendy relations 

The logR e -<yu> e diagrams for the galaxies of A2163B are 
shown in Fig. [TO] for the R and K bands. We note that the I- 
band data provide almost the same spectral information as the 
R-band photometry, and therefore are not discussed in the fol- 
lowing. For each band, we excluded the disks (n < 2) and the 
galaxies for which the total magnitude estimated by the struc- 
tural parameters is fainter than the completeness magnitude. 
This selection results in N = 52 and N = 53 galaxies in the R 
and K bands, respectively. For what concerns the NIR data, we 
also excluded two galaxies that have small radii and low sur- 
face brightness and for which, therefore, structural parameters 
are possibly affected by a larger uncertainty. The K-band sam- 
ple consists of N = 5 1 objects. In order to describe the log R e - 
<fi > e sequence, we consider the following equation for each 
waveband M: 

om+jSm- log (5) 



where R e is the effective radius in kpc, au and /3m are the zero- 
point and the slope of the relation. The evolution of the slope, 
of the zeropoint and of the intrinsic dispersion of the KR has 
been recently studied in the optical by LBM03b, who showed 
that these properties do not change significantly over redshift 
at least up to z ~ 0.6. They found (3 = 2.92 + 0.08, an intrin- 
sic dispersion of 0.4 + 0.03 mag/arcsec 2 , and a KR zeropoint 
of 18.95 + 0.08 mag/arcsec 2 in the R band at z ~ 0.21. In or- 
der to derive the values of a and /3 for A2163B, we applied 
the BMLS fitting procedure (see LBM03b), which allows the 
selection criteria to be taken into account, providing un-biased 
estimates of the KR coefficients. The fitting results are shown 
in Fig. [TO] For what concerns the dispersion of the KRs, its 
value turns out to be consistent between the R and K bands, 
amounting to 0.65 + 0. 12 in <jj > e and to 0.22 + 0.04 in log R e . 
The R-band slope of the logR e -<yu> e relation is fully con- 
sistent with that of LBM03b, in agreement with their finding 
that this coefficient does not evolve significantly over redshift. 
By adopting /3 = 2.92 in the fit, we obtain a R = 19.16 + 0.1 
and, therefore, also the R-band zeropoint turns out to be in full 
agreement with the result of LBM03b, once that measurement 
errors are taken into account. By subtracting in quadrature from 
the observed scatter around the KR the amount of dispersion 
due to the measurement uncertainties 13 on logR e and <fi> e , 
we obtain the following estimate of the intrinsic dispersion of 
the KR: 0.6 + 0.13 mag/arcsec 2 (0.21 + 0.04 dex in R e ). This 



To this aim, we took into account the correlation of the uncertain- 
ties on log R c and </i > c . 
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Fig. 9. Distributions of V - R and V - K colour gradients. The mean values are marked by a solid line. 



value is a little larger than, but consistent with the estimate of R e and </x> e due to internal colour gradients of galaxies. By 
LBM03b. using Eqs. l,|4]and[5l we obtain the following relation: 

<fi >f = <fi >f -5 ■ (log Rf - log Rf) - (R T - K T ) - 
(1 + A)" 1 • [a R -T + 2.5 ■ A ■ log In - (5 - f3 R ) ■ 

■logj§+ (^ R + 5A)-logR e K ] = 



(6) 



T , , „ . . , ■ , , .„„ where log ^ is the mean logarithmic ratio of R- and K-band 

In order to compare the coefficients of the optical and NIR R e 

logR e -<jU> e sequences, we have to take into account (1) the effective radii, T and A are the zeropoint and the slope 14 

variation of the galaxy colour along the sequence due to the 14 We note that Eq.|6lwas derived assuming that eq. 1 still holds 

colour magnitude relation and (2) the waveband dependence of once that aperture colours are replaced by the total galaxy colours, 
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Fig. 10. Kormendy relations in R and K bands. Dashed lines mark the completeness magnitudes and the surface brightness cut 
(right panel) of each sample. The typical uncertainties on the effective parameters are indicated by the ellipses in the lower left 
corner of each panel. The ellipses mark the contours that enclose a probability of 68% for a 2D normal deviate. 
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Fig. 11. Ratios of optical to NIR effective radii and colour gradients are plotted versus the logarithm of the K-band effective 
radius. The dashed lines mark the mean values of log R^/Rf and grad(V - K) in the left and right panels, respectively. Data were 
binned in order to have the same number of points in each bin. The solid lines mark the best-fits to the data, while the dot-dashed 
lines indicate the lcr confidence intervals. 



of the R - K vs. Kj colour magnitude relation, a' K 

(1 +A)- 1 ■ [or - T + 2.5 ■ A • log 2n - (5 -00 • log(R*/R e K ) 

p K = (1 + A)" 1 • (f3 R + 5A), and y' K = - (5 -fig) • (1 + A)" 1 . 
We note that the value of a' K corresponds to the zeropoint of 
the NIR KR, while fi' K is equal to /?k only if the third term of 
the last member of Eq.|6] vanishes, that is if the ratio of optical 
to NIR effective radii and, therefore, the internal colour gradi- 
ent of galaxies does not vary systematically along the log R e - 
<fi > e sequence. By using the coefficients of the CM relation 
derived in Sec. 3.1, the mean ratio of optical to NIR effective 
radii (~ 0.16), and the R-band KR coefficients (a = 19.16 and 
P = 2.92), we obtain a' K = 15.57,/^ = 2.76 and y' K = -2.23. 
The value of a' K is in good agreement with that of qtk, while 



the value of f3' K is only marginally consistent with /3k, with a 
difference of 2.3cr. Since the value of y' K is negative, this dif- 
ference could be explained if the R^/R^ ratio would become 
smaller for larger values of R^, implying that larger galaxies 
have a less steep colour gradient. Fig. II 1 blots log(R^/R^) and 
grad(R - K) as a function of log R^. Although the figure seems 
to favor the trend suggested by the difference of f3' K and [3k, the 
error bars do not allow such a trend to be definitively proven. 
We point out, however, that the previous results are in clear dis- 
agreement with a possible steepening of colour gradients with 
galaxy size, since it would make the discrepancy between /3k 
and P' K strongly significant. This point will be further addressed 
in Sec. |9] 



Rt - Kt. We verified, in fact, that the values of T and A do not change 
significantly by using Rron magnitudes to estimate the galaxy colours. 
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Fig. 12. Distribution of galaxies in the space of structural parameters. Upper-left: logR e -logn relation. The ellipses is in the 
upper plots are the 68% probability contours. The solid and dot-dashed lines are the best-fit and the \<j confidence lines, 
respectively. The region affected by the magnitude selection is filled by a dashed pattern. Upper-right: edge-on projection of the 
PHP with the best-fit relation marked by a solid line. Lower-left: 'long' edge-on projection of the PHP. Lower-right: face-on 
view of the PHP. The shaded region is delimited on the right by the line that results from the intersection of the plane K T = const, 
with a plane parallel to the PHP but with zeropoint c + 2 ■ <x c . 



8. The NIR Photometric Plane at z ~ 0.2 

Fig-Elplots different projections of the distribution of galax- 
ies in the space of the NIR structural parameters 15 . As shown in 
Figll2h. effec tive radii are correlated with Sersic indices (e.g. 
lYoung and Curriei200lh . in the sense that galaxies with a larger 
size tend to have a surface brightness profile more peaked in the 
center. By looking at Fig.^3 we see that, taking into account 
the dispersion of the KR, the completeness magnitude line be- 
gins to intercept the log R e -<yU> e relation at logR e ~ 0.3, im- 
plying that galaxies in the shaded area of Fig I12h start to be 
lost from the present sample because of the magnitude selec- 
tion. We applied a bi-weight least square fit to the data with 
log R e > 0.3, minimizing the rms of the residuals with respect 
to logn. The fit gives logn = (0.31 +0.05) -log R e + (0.56 + 0.1), 
with an rms that amounts to ~ 0.15dex (~ 35%) in n and 
~ 0.5dex in logR e , which is much larger with respect to that 
obtained for the KR (~ 0.22 + 0.04 dex). We note that the log n 
regression is quite insensitive to selection cuts with respect to 
log R e , and, in fact, the above coefficients do not vary signifi- 
cantly including all the galaxies in the fit. 



15 The waveband dependence of the Photometric Plane will be inves- 
tigated in a forthcoming paper, by using a larger sample of galaxies. 



The other panels of Fig.lT2"lshow different projections of the so- 
called Photometric Plane (PHP). For comparison with GRA02, 
we adopted the following representation of the PHP: 

log R e = a ■ log n + b <\i > e +c, (7) 

which is analogous to the usual equation of the Fundamental 
Plane (FP), once the velocity dispersion term is replaced with 
logn. Following GRA02, the coefficients a and b were ob- 
tained by a least square fit minimizing the rms of the residu- 
als to logR e . We also took into account the correlation of the 
uncertainti es on the structural parameters by using the MIST 
algorithm (La Barbera. Busarello and Cap accioli 12000). The 
fitting results are a = 0.6 + 0.13, b = 0.235 + 0.02 and 
c = -4.5 + 0.25, with a residual rms of ~ 0.15 dex in Re. In 
order to correct these values for the magnitude cut, we general- 
ized to the 3D case the procedure adopted in the previous sec- 
tion for the KR fit (see LBM03b for details). We found that the 
bias in the b coefficient is negligible, while it amounts to +30% 
and -4% for the coefficients a and c, respectively, and to +13% 
for the dispersion in the log R e variable. After correction, the 
fitting coefficients become a = 0.8 + 0.2, b = 0.235 + 0.02 
and c = -4.70 + 0.026, with a corrected logR e scatter of 
~ 0.17 dex, which is smaller by ~ 14% than that of the KR. The 
'corrected' values of a, b and c were used to obtain the edge-on 
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sections (Fig. 112b. c) and the face-on projection (Fig.ll2H) of 
the PHP. 

The slopes and the scatter of the NIR PHP at z ~ 0.2 turn out 
to be in remarkable agreement with those derived by GRA02 
for galaxies in the Fornax and Virgo clusters: a = 0.86 + 0.13, 
b = 0.228 ± 0.036 and o- logRe ~ 0.17 dex. We point out that 
the magnitude selection significantly affects the logn coeffi- 
cient and that, therefore, it would be very interesting to use a 
sample of galaxies spanning a larger magnitude range in order 
to better constrain this coefficient. The selection effects can be 
further illustrated by looking at the face-on projection of the 
PHP in Fig. I12H . Due to the finite dispersion of the galaxies 
around the PHP, its intersection with the plane Kj = cost re- 
sults in a 'strip', whose 'upper' border marks the limits of the 
dashed area in Fig. I12H . As consequence, the distribution of 
galaxies in this area is not 'complete'. 

In order to estimate the intrinsic dispersions of the log R e -log n 
relation and of the PHP, we subtracted to the observed disper- 
sion in log R e the amount of scatter expected from the measure- 
ment errors on the observed parameters, taking into account the 
covariance terms between the uncertainties on logR e , </i> e 
and log n. The intrinsic dispersions of the log R e -log n relation 
and of the PHP turn out to be ~ 0.45 dex and ~ 0.16 dex, 
respectively. Interestingly, we find that the measurement er- 
rors account only for few percents of the observed dispersion 
around the PHP, because, as shown in Fig. 112b. their correla- 
tion is almost parallel to the plane, in the same way as for the 
ICR. 

9. Discussion and conclusions. 

We have studied the galaxy population in a 5' x 5' field around 
the X-ray extended source A2163B, at about 6.5' North the 
cluster of galaxies A 2163 (z = 0.201). The BVRIK catalogue 
is presented, including total magnitudes, colours, photometric 
redshifts and the RIK structural parameters (effective radius, 
mean surface brightness and Sersic index). 

The K-band luminosity density map of A2163B presents 
a main central overdensity of radius ~ 0.2 Mpc, and a sec- 
ondary structure ~ 0.4 Mpc from the center. The colour mag- 
nitude diagrams show a sharp red sequence whose slope flat- 
tens at longer wavelengths: galaxy colours become bluer by 
~ 0.1 mag per magnitude in the B - K vs. Kt diagrams, and 
by ~ 0.04 mag into the I - K vs. Kt plane. These values agree 
with the predictions of a mass-metallicity model of the CM 
relation (KAB98). By introducing a suitable procedure to ac- 
count for the Milky Way reddening, we find that the zeropoints 
of the CM sequence are consistent with those expected for a red 
galaxy population at redshift z c ~ 0.2, implying that A2163B 
is a cluster of galaxies likely involved in a merging event with 
A 2163. The value of z c is also confirmed by the analysis of the 
K-band LF and by the photometric redshift technique, which 
is used to select cluster members. The luminosity segregation 
of A2163B is investigated by means of the K-band LF. The 
global LF is well described by a single Schecter function with 
a faint end slope of a = -1.2, in agreement with other studies 
of the NIR LF both in clusters and in the field (e.g. lAndreonl 
2001, Kochanek et al. 2001). On the other hand, the faint end 



slope depends on environment, varying from about -1 in the 
two higher density clumps, to about -1.4 into the outer region, 
implying that the population of dwarf galaxies is less abundant 
into the denser environment. This behaviour of the LF suggests 
that some mechanism, such as tidal disruption or cannibalism, 
is working in the higher density environments, decreasing the 
number of faint galaxies (Haines et al. 2003, submitted). 

Surface photometry has been derived in the R, I and K 
bands for N ~ 60 cluster galaxies, and has been used: 

i) to compare galaxy parameters at optical/NIR wavelengths, 

and to estimate the optical and NIR internal colour gradi- 
ents; 

ii) to perform a comparative optical/NIR analysis of the KR; 

iii) to derive the PHP of spheroids at z ~ 0.2 in the K band. 

The characteristics of the light distribution in galaxies are the 
same between the optical wavebands, while a significant dif- 
ference is found when comparing optical to NIR parameters. 
Effective radii become a factor 1.45 ± 0.15 greater from the 
K to the R band, implying that the light profile in galaxies is 
much more concentrated in the NIR than in the optical. On 
the other hand, surface brightnesses show a waveband depen- 
dence that is fully explained by the behaviour of R e and by 
the R - K galaxy colours. The value of R^^/R^ 111 is con- 
sistent with what we found in previous works on the opti- 
cal/NIR structural parameters of cluster galaxies at z ~ 0.3 and 
~ z ~ 0.6 (see LBM02 and LBM03a) , and is within the r ange 
of results repor ted in previous works. Nels on et alJ (l2002h and 
iRembold et ail d2002h found R^/R^ ~ 2 by analysing 
brightest cluster galaxies at 0.4 < z < 0.8 and nearby early- 
types, re spectively, while a lower value, RP^/R^™ ~ 1.2, was 
found bv lPahre. de Carvalho and Diorgo vski ( 1998). 
The population of spheroids in A2163B has null colour gra- 
dients at optical wavelengths (V - R restframe) and optical- 
NIR gradients that are significantly negative (~ -0.48 + 
0.06 mag/dex), consistently with the quoted differences among 
the effective radii. As shown by LBM03a, the optical-NIR 
gradients do not evolve significantly with redshift, amount- 
ing to -0.4 + 0.1. As detailed in that work, the redshift de- 
pendence of both optical-NIR and UV-optical colour gradi- 
ents strongly constrains the age and metallicity of the inner 
and outer galaxy stellar populations. Using the colour gradi- 
ent model of LBM03a (see model Zl, Tab. 6), the quoted 
colour gradients imply a metallicity gradient in spheroids of 
-0.2 + O.ldex per decade of radius. This result is in agree- 
ment with that reported from other studies of colour gradients 
for both nearby and intermediate redshi ft galaxies (Sagli a~t alJ 
l2000Hldiart.'M ichard a nd Pachecol2002l) . 

Monolithic collapse models predict that early-type galaxies 
are more metal rich in the center than in the outskirts. During 
galaxy formation, the gas dissipates its kinetic energy carry- 
ing the heavy elemen ts ejec ted from the evolved stars into the 
inner region dLarsori Il974l) . Numerical simul ations predict a 
meta l abundance g radient of about -0.5 (e.g. Carl berg|ll984[ 
and lKawatalboOlt see their models B1/B2 which correspond 
to the magnitude range of our samples). This value is not con- 
sistent with our results, although some caution is needed due 
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to scatter present both in the ob served data and in the simu- 
lation results (see Kawata 2001). The difference of observed 
and predicted metallicity gradients suggests that gas dissipa- 
tion has a minor role in monolithic collapse models. We note, 
however, that dissipation is required in order to reproduce the 
mass-metallicity relation of early-type galaxies and their mod- 
erate orbital anisotropies. We argue, therefore, that some tun- 
ing of the parameters driving the physics of gas dissipation 
would be required in order to fit all the observed properties into 
the monolithic scenario. A likely explanation of the observed 
metallicity gradients resides in the merging processes, which 
mix the SPs inside galaxies. It has been sh own that both dis- 
sipationless merging of spheroidal systems JWhiteLll980l) and 
dissipative merging of disk galaxies (Bekki & Shiova, 2001) 
produce more shallow stellar population gradients. 

The slope and the zeropoint of the optical log R e -< fi > e re- 
lation are consistent with those obtained in LBM03b. For what 
concerns the NIR KR, we compared its slope and zeropoint 
with the optical values by taking into account (i) the colour 
magnitude relation and (ii) the mean value of the ratio of op- 
tical/NIR effective radii. These factors fully explain the differ- 
ence among the optical and NIR zeropoints. The most inter- 
esting result is that the slope of the R-band KR, corrected for 
the effects (i) and (ii), is lower by ~ 2.3<x than the K-band 
value. We show that this difference can be explained if the ratio 
of optical to NIR effective radii (the optical/NIR colour gra- 
dient) decreases (flattens) for larger galaxies. A direct inspec- 
tion of the log R^/Rf and the grad(V - K) vs. log R^ diagrams 
shows that such trend exists. Our data do not favour a steep- 
ening of colour gradients with galaxy size, that would imply 
a larger discrepancy between the R- and K-band KR slopes. 
A steepening of colour gradients with galaxy luminosity (size) 
is a natural expectation of the monolithic collapse model, due 
to the galactic wind mechanism JLarsori Il974l) . The galactic 
wind blows early in less massive galaxies, preventing gas dis- 
sipation to carry heavy metals in the center, and producing, 
therefore, a less steep gradient in these systems. A steepening 
of colour gradients with galaxy luminosit y has been unfruit- 
fully looked for by various studies (e.g. [P eletier et al. 1990, 
iBartholomew et al. lEoOll) . Recently. iTamura and OhtTJ2003l) 
measured the optical colour gradients for N = 36 galaxies at 
z ~ 0. They found that the optical colour gradient of very bright 
early-types (L > L*) can steepen with galaxy luminosity, while 
an opposite trend was suggested f or fainter galaxies. A similar 
bimodal behaviour was found by (Bal cells and Peletien. 11994) 
for the colour gradients of early- type spiral bulges. We note 
that our data do not allow us to distinguish the behaviour of 
colour gradients for the very bright early-types, while they al- 
low a larger magnitude range in both optical and NIR to be 
analyzed. In hierarchical merging models, more massive galax- 
ies form from larger disk systems, which have stronger inter- 
nal gra dients. Since merging dilutes the stellar population gra- 
dients JWhiteL fT980). if gala xies with larger size have ex peri- 
enced a higher merging rate ( Cap accioli. Caon and D'Onofriol 
119921) . their colour profile can be as flat as or even less steep 
than those of smaller systems. The data of A2163B seem to 
favour this interpretation. 



Finally, we have analyzed the so-called Photometric Plane 
(PHP) of spheroids at z ~ 0.2. We find that galaxies follow a 
bivariate relation between K-band effective parameters and the 
logarithm of the Sersic index: log R^ oc a • log n K + b- <[i>f. 
Accounting for the correlations among the uncertainties on 
structural parameters as well as for the selection effects, we 
find a = 0.8 + 0.2 and 0.235 + 0.02, and an observed disper- 
sion of ~ 0.17 dex in logR e . These values are consistent with 
the optical PHP coefficients found by GRA02 for early-types 
in the Fornax and Coma clusters. We have discussed how the 
magnitude cut affects the distribution of galaxies in the space 
of the structural parameters: the logn term of the PHP turns 
out to be significantly affected by the magnitude selection, in- 
dicating that it would be very interesting to analyze a larger, 
deeper sample of galaxies to further constrain the coefficients 
of this relation. Interestingly, we find that the correlation of the 
uncertainties on the structural parameters is almost parallel to 
the PHP, analogously to what happens for the KR. This implies 
that the intrinsic scatter of the PHP amounts to about 0.16 dex, 
which is only few percents smaller than the observed scatter, 
but ~ 13% lower than that of the KR (~ 0.21 dex). Therefore, 
although the observed scatter around the plane seems to be sim- 
ilar to that of the spectroscopic Fundamental Plane (~ 0. 15 dex) 
as found by GRA02, the present data indicate that its intrinsic 
dispersion can be significantly larger. This point will have to 
be further analyzed by studying samples of galaxies for which 
both velocity dispersions and Sersic indices are available. 

Appendix A: Data reduction 

A.1. BVRIdata 

Both the images of the first and of the second run were reduced 
by following the same procedure. After bias subtraction, 
the images were divided by a flat-field frame obtained by 
combining twilight sky exposures. Since this procedure did 
not fully remove low frequency variations of the background 
across the images (at level of ~ 5%), we obtained a super-flat 
correction for each band by median combining the correspond- 
ing scientific exposures. A polynomial fit was applied to each 
super-flat frame and the fitted surfaces were used to correct 
the scientific exposures. This procedure reduced background 
variations in the final images at the level of ~ 0.5%. After 
flat-field, the cosmic rays in each exposure were detected by 
using the IRAF task COSMICRAYS, and a mask frame was 
constructed including hot pixels and the bad columns. The 
images of each run were then combined by using the IRAF 
task IMCOMBINE. Since only Run II was photometric, the 
images of Run I were suitably scaled and combined with those 
of Run II. The scaling factors were computed from the mean 
difference of the magnitudes of bright unsaturated stars in the 
observed field. The seeing FWHM is 1.5", 1.2", 1.0" and 1.0" 
for the B, V, R, and I bands, respectively. The images were 
calibrated into the Johnson- Kron-Cousins p hotometric system 
by using standard stars from Landolt ( 1992). Since comparison 
fields for different airmasses were not available, we adopted 
the extinction coefficients typical for the site of La Silla: 
A B = 0.21, A v = 0.13, A R = 0.09, Ai = 0.05 mag/airmass. 



F. La Barbera et al.: The cluster of galaxies A 2163B. 



17 



Instrumental magnitudes for the Landolt stars were computed 
with in an aperture of d iameter 10" by means of S -Extractor 
jBertin & Arnoutt| 19961) . For each band, in order to obtain the 
zeropoint and the colour term of the photometric calibration, 
we performed a robust least square fit. Colour terms were 
found to be small and were neglected, while the zeropoints 
turned out to be ZP B = 24.31 + 0.01, ZP V = 24.92 + 0.01, 
ZP R = 25.06 + 0.01, ZP : = 24.41 ± 0.01 mag (scaled to Is 
exposure time). 

A.2. K-band data 

For each dithering sequence, the exposures were dark sub- 
tracted and corrected for flat-field, using a super flat frame ob- 
tained by median combining all the images taken during the 
night. After this procedure, the magnitudes of the standard stars 
showed a rms variation across the chip of ~ 0.04 mag. To 
achieve a better accuracy, we obtained an illumination correc- 
tion frame by measuring the magnitude of a standard star at 
different positions in a 5 x 5 grid on the frame. All the im- 
ages were divided by the illumination correction frame, which 
allowed the low frequency component of the flat-field to be cor- 
rected at better than 1%. Since sky subtraction is a very trou- 
blesome step for the reduction of NIR data, particularly in high 
density regions mostly populated by bright early-type galax- 
ies with extended halos, this point was carefully dealt with by 
a two step procedure. First, each exposure was sky subtracted 
by computing the sky frame from the median of the six clos- 
est frames along the sequence. The images were then regis- 
tered with integer shifts and combined, with a sigma clipping 
algorithm for cosmic ray rejection. This procedure, however, 
overestimate the sky level in the extended halos of galaxies, 
and therefore, it is not suitable to derive the surface photom- 
etry. For this reason, we reiterated the sky subtraction, pro- 
ceeding as follows. The first step images were used to ob- 
tain mask frames for the sources in the field, by running S- 
Extractor with the checkimage OBJECTS option. For each se- 
quence, the mask was expanded in order to 'cover' the galaxy 
halos and was de-registered to each dithered exposure, by in- 
cluding also hot and bad pixels. The sky frames were esti- 
mated by an average of the six closest frame for each se- 
quence exposure, rejecting masked pixels. The images were 
then sky subtracted and combined with IMCOMBINE 16 (using 
the SIGCLIP option), resulting in a final image with a seeing 
FWHM of ~ 0.9". The photomet ric calibration w as performed 
into the Ks standard filter JPersson et all deriving the 

instrumental magnitudes of the standard stars within an aper- 
ture of diameter 8". The airmass correction was performed by 
using an extinction coefficient Ak = 0.06 mag /airmass, which 
was derived by comparing the magnitudes of bright objects in 
the field for different airmasses. The K-band zeropoint turned 
out to be ZPk = 22.422 + 0.015 mag (scaled to 1 s exposure 
time). 



16 The images taken under non-photometric conditions were suitably 
scaled to the others. 



Appendix B: Deriving z c from the CM relations 

The zeropoints Tk, of the CM relations (Eq.Q were fitted with 
the colours calculated for an old, passively evolving SP, as a 
function of the redshift z c and of the reddening E(B - V). To 
account for the metallicity change along th e CM sequence , we 
proceeded as follows. As shown by Merlu zzi et all (|2003), the 
CM relation is reproduced by the metallicity-luminosity rela- 
tion: 

log(Z/Z ) = (-0.097 + 0.005) • M v + (-2.09 + 0.09). (B.l) 

where My is the absolute V-band magnitude. Using the slope 
of the CM relation for nearby early-types from Bower et 
al. (1992), the above relation allows the absolute K-band mag- 
nitude to be estimated for a given metallicity Z. For Z = Z , 
we obtain K(Z ) = K = -24.74 + 0.15. In order to pre- 
dict the value of K(Z ) as a function of redshift (K (z)), we 
have to take into account the luminosity distance term, and the 
K+E corrections. To obtain these c orrections, we used galax y 
templates from the GISSEL00 code feruzual & CharloJl993l) . 
with a Scalo IMF, an exponential star formation rate exp(-t/r) 
and solar metallicity. These models were also used to predict 
(1) the optical-NIR galaxy colours as a function of the redshift 
z c , and (2) the absorption coefficients An and Ak as a func- 
tion of z c and E(B - V), using the extinction law from SFD98. 
The values of z c and E(B - V) were derived by minimizing the 
following expression: 

* 2 [z c ,z f ,E(B - V)] = J] [r Ko(z ) - (N - K) t + (A N - A K )] 2 (B.2) 

M 

where (N - K) t are the template colours, and Zf is the forma- 
tion redshift, which was allowed to vary in the range [1.2, 10]. 
The uncertainties on z c and E(B - V) were obtained by shifting 
the values of Tk^z) according to the relative uncertainties and 
re-iterating the y 2 minimization. We note that the uncertainties 
on Tk f also include the contribution of the errors in photometric 
zeropoints, which were added in quadrature to the uncertainties 
on galaxy magnitudes (see Sec. 12. 2\ . We found that the esti- 
mates of z c and E(B - V) do not change significantly by adopt- 
ing either t = 0.01 Gyr, (Simple Stellar Population models), or 
t = 1 Gyr (models with a more protracted star formation rate). 
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